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Magnetic field plays a very important role in many astronomical phenomena at various scales 
of tlie universe. It is no exception in tlie early universe. Since the energy density, pressure, and 
tension of the primordial magnetic field affect gravitational collapses of plasma, the formation of 
- - - seeds for large scale structures should be influenced by them. Here we numerically investigate the 

effects of stochastic primordial magnetic fleld on the seeds of large scale structures in the universe 
. in detail. We found that the amplitude ratio between the density spectra with and without PMF 

■ {\P{k)/Po{k)\ at fc > 0.2 Mpc~^) lies between 75% and 130% at present for the range of PMF 

' ^ I strengths 0.5 uG < B\ < 1.0 nG, depending on the spectral index of PMF and the correlation 

O ■ between the matter density and the PMF distributions. 
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INTRODUCTION 



, The possible existence of a primordial magnetic field (PMF) is an important consideration in modern cosmology. 

The origin of PMFs has been intensively studied [11, [3, H, 0, M, H S 3 ■ Although it seems unlikely to survive an epoch 
I of inflation, it is conceivable that large-scale magnetic fields and magnetic inhomogeneities could be generated at the 
■ end of that era or in subsequent phase transitions @. Studies of magnetogenesis are partly motivated by the need 
" to explain the origin of large-scale magnetic fields which are observed in galaxies or in clusters of galaxies [l^]- If 
' these ma gne tic fields have the primordial origin, the PMF should have influenced a variety of phe nomena in the early 
Q-i universe [111, [11 III [11, 111, III [O, III] ■ Several semi-analytic and numerical studies [H 111 111 IH, 113, [Hi pointed 
Q I out that the effect of the PMF is one of the new physical processes in the early universe. 

5-H ' Large-scale magnetic fields are now considered to be standard components. Even at high redshift, the existence 
c/3 , of dynamically significant magnetic fields is suggested from observations of Faraday rotation associated with high 
' redshift Lyman-a absorption systems [lo| . However their picture in the early universe is not clear and is still a matter 
^ ' of debate. If dynamically significant large-scale magnetic fields were present in the early universe, they must have 
affected the formation and evolution of the observed structure, and some signatures of their past should be included 
in their structure and spectrum. Indeed, high resolution Faraday rotation maps and the study of extragalactic cosmic 



. rays provide direct observational support for this point of view [21| . 
" " ' In the evolution of baryonic and dark matter at large scales, accretion shocks form in the infalling flows towards 
the growing nonlinear structures such as sheets, filaments and clusters (e.g., [H H^). In fact, some evidence of 
these accretion shocks has been detected in radio relic sources [IJ] . The properties of the shocks depend on the power 
spectrum of the initial perturbations of baryonic and dark matter on a given scale as well as the background expansion 
in a given cosmological model. Therefore, if there existed dynamically significant large-scale magnetic fields in the 
early universe, the dynamical influence of magnetic fields on the spectrum of perturbations and thus on flow motions 
should not be ignored. From this point of view, investigations of the physical process in the early universe with PMF 
provide important suggestions for studies of the formation and evolution of large scale structures (LSS). 

In this paper, we assume the existence of a large-scale PMF and analyze the magnetic effects on density in- 
homogeneities numerically. Effects of PMF on density fields, especially on cosmic microwave background (CMB) 
anisotropics at the photons' last scattering surface (PLSS), were studied by considering the scalar-type component of 
energy momentum tensor in PMF , or including a new analytical magnetic power spectrum source due to a Lorenz 
force without previous approximations. Here we include both effects consistently and extend their previous studies 
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by expanding the analysis from the epoch of creation of CMB anisotropics toward the present epoch. In particular, 
we consider the different correlation between the PMF and the matther power specturm, and investigate their effects 
on the density fields at large scales in the presence of PMF. 

Throughout this paper we fix the cosmological parameters as follows [27|: h — 0.71, fib — 0.044, ^cdm — 0.226, 
Us — 0.93, and Tc — 0.10 in flat Universe models (thus Q.\ — 1 — VLi, — flcDM — 0.73), where h denotes Hubble 
parameter in units of 100 km/s/Mpc, fit and Qm are the baryon and cold dark matter densities in critical density 
units. Us is the spectral index of the primordial scalar fluctuation, and Tc is the optical depth of Compton scattering. 



II. COSMOLOGICAL PERTURBATIONS WITH STOCHASTIC PRIMORDIAL MAGNETIC FIELD 



A. Cosmological MHD 

Let us consider the PMF created by some effects during the radiation-dominated epoch. The energy density of 
the magnetic field is treated as a first order perturbation in a flat Friedmann-Robertson- Walker (FRW) background 
cosmology. The electromagnetic tensor has the usual form 
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where Ei and Bi are the electric and magnetic fields. The energy momentum tensor for electromagnetism is 

^"^[EM] = ^ - Id'^^F^sF'''^ . (2) 

The Maxwell stress tensor, consisting of the space-space components of the energy momentum tensor of the electro- 
magnetic field, is 

Tik ik 
[EM] = cr = 

J_J_/^«^fc + B^Sfe_i^»fe(£;2 + s2)\ (3) 

Oj T:7r I ^ J 

Within the linear approximation, the magnetic field evolves as a stiff source. Therefore we can discard all back 
reactions from the magneto hydrodynamic (MHD) fluid onto the field itself. The conductivity of the primordial 
plasma is very large, and the field is "frozen- in" [1J| . This is a very good approximation for the epochs in which we 
are interested. Furthermore, we can neglect the electric field, E ^ 0, and also decouple the time evolution of the 
magnetic field from its spatial dependence, i.e., B{t, x) = BQ{x)/a? for very large scales. In this way we obtain the 
following expressions, 

B^ 

r'^°[EM] 1 F , (4) 

T'°[EM] 7'°*'' [EM] = , (5) 

-T^'[EM] = 'T*' - ^(25^5" - ^''^')- (6) 
First and second terms in Eq.® arc magnetic tension and pressure, respectively. 



B. Evolution Equations of Cosmological perturbations 

Combining the Einstein equations and linearizing them, we obtain the perturbation evolution equations. In order 
to consider the effect of the PMF, we should add the electromagnetic energy momentum tensor to that of standard 
cosmic fluids, 

r"/' = r"^[FL]+r"^[EM]- (7) 

We assume that the baryonic matter is representable as a perfect fluid and neglect the anisotropic pressure pertur- 
bations. As an initial condition, we consider adiabatic perturbations and neglect entropy perturbations initially. The 
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line element in the conformal synchronous gauge for a flat Friedmann-Robertson- Walker (FRW) background is given 

by 

ds^ = a'^{T)[-dT^ + {d^j + h^j)dx'dx^], (8) 

where is the spatial coordinate, a(r) is scale factor, r is the conformal time defined by dr — dt/a^r), hij are metric 
perturbations around the background spacetime, and the speed of light is set to unity. We will be working in the 
Fourier space in this paper. We introduce two fields h{k, r) and ri{k, r) in fc-space and write the scalar mode of hij as 
a Fourier integral 



hij{x,T) = J 



kikjh{k,T) + 6?7(fc, r) ( kikj ~ -6, 



1 



(9) 



where k is the wave number in the Fourier space and k is kk and fe is a unit vector of wave number. The linearized 
perturbation equations are obtained from the Einstein equations up to first order [25L Isil [33 | : 

k^T] - ^Hh = AnGa^ST", (10) 

fc^r) = AnGa^ik^T^, (11) 

h + 2Hh-2k^rj = SnGa^ST^, (12) 

h + 6ii + 2H{h + 6ri) -2k^r] - ''''-^"^ I '^x. \lxrri J ^ 
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= 247rGa2{Z[EM:S](fc)-(p + p)^} (13) 

where 

- 5Tq = 6p = (5/?[FL] + ^^[EM] , (14) 

ik^6T° = {p + p)v , (15) 



st; = <5r;pL] + ^t;^^ , (i6) 

Z[EM:S]{f^) is a scalar part of the magnetic shear stress, and {p + p)(j is a fluid shear stress. We assume that the PMF 
Bq is statistically homogeneous, isotropic and random. For such a magnetic field, the power spectrum can be taken as 
a power-law S{k) =< B{k)B*{k) >oc fc" [14] where n is the power-law spectral index of the PMF. The index n can be 



either negative or positive depending on the physical processes of the creation. From ref. |14| , a two-point correlation 
function for PMF is defined by 

B\k)B^\U)) = ^—l^-^k-P^^{k)S{k-k'), k<kc, (17) 



where 



P'\k) = 5^'-^. (18) 

B\ is the magnetic comoving mean-field amplitude obtained by smoothing over a Gaussian sphere of comoving radius 
A, and k\ = 27r/A (A = 1 Mpc in this paper). The cutoff wave number kc in the magnetic power spectrum is defined 
by 0, 

4^(p+p) Jo '^^ T' < '^dec (^9) 
fcc "(■^doc), T > Tdoc, 

where 1^ is the mean free path of photons, and Tdcc is the time of the decoupling of photons from baryons(see appendix 
B). Evaluating the two-point correlation function of the electromagnetic stress-energy tensor in the Fourier space, we 
obtain the power spectrum of PMF energy density, Lorenz force and shear stress as the following (see appendix) 

|£^[EM:S](fc,T)p5(fc- fc') = (r(fc,T)[EM:Sl]r*(fc',r)[EM:Sl]) , (20) 
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|n[EM:S](fc,r)|2<5(fc-fc') = ^^((r(fc,T)[EM:Sl] -T(fer)[EM:S2]) (r*(fc',T)[EM:Sl] -T*(fc',r)[EM:S2])), (21) 



and 



\Z{k)[EM:S]\^S{k^k') = ^ ^0r(fc,r)[EM:Sl] -r(fc,r)[EM:S2]) t)[eM:S1] " ^)[EM:S2]) ^ , (22) 

respectively. Explicit expressions for the ensamble averages to evaluate the above spectra, in the case of power law 
stochastic magnetic field, are given as (also see appendix) 
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(%M:Sl](fc)?'[*EM:S2](^)) + (%M:S2](fc)r[EM:Sl](^)) 

2 



1 J (27r)"+8 B2 



(2^)7a8 I 2fc5:+3 r (s±^) 
3 



c?fc fc 



IUn+3 



{kk'fn{n + i) 



{{k + k'Y+^ -\k-k'y'+'^} 



+ 



fc2fc'3n(n + 2)(n + 4) 
3 

fc3fc'4n(n + 2)(n + 4)(n + 6) 



(r[EM:S2] (fc)r[EM:S2] (^)) 



||fc-/c'|"+4 + (fc + fc')»+4} 



k^k'^n{n + 2){n + 4) 

{(fc + fc')"+''- i/c-fe'r^*^} 



{(fc + fc')"+'*-l^-fcT+''} 



(24) 



1 j (27r)"+« ^2 



(2^)7a8 2fc 



n+3 p /' n+3 



dk'k' 



{kk'){n + 6) 



(fcfc')%(n + 2)(n + 4) 
3 



(fcfc')2(n + 6)(n + 8) 



{(A: + fc')"+^-|A:-A:T+*} 



{(fc + fc')"+''-|'t-fc'l"+'} 



(25) 



The evolutions of fluid variables can be obtained by imposing the conservation of energy-momentum, which is a 
consequence of the Einstein equations 



(26) 



This leads the following equations in /c-space: 
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(27) 

(28) 



where w = p/ p. In the continuity and Euler equations for the scalar mode, we can just add the energy density and 
pressure of the PMF to the energy density and pressure of cosmic fluids, respectively. Since baryon fluid behaves as a 
nonrelativistic fluid in the epoch of interest, we may neglect w and 6pb/ Sph, except the acoustic term Cgk^Sb- Also, the 
shear stress of baryons is far smaller, and we can neglect it [1^. Since we concentrate on scalar type perturbations in 
this paper, we do not consider the magneto-rotational instability from the shear stress of the PMF and baryon fluid 
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From equations ((27)) and p8|) . and by considering Compton interaction between baryons a nd p hotons, we obtain 
the same form of the cvohition equations of photons and baryons as in previous works |26l . [sil . |32| |. 
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(33) 



where Tig is the free electron density, ctt is the Thomson scattering cross section, and a-^ of the second term on the 
right hand side of equation ((3T|) is the shear stress of the photon with the PMF. Since n < is favored by constraints 
from the gravitational wave background [38] and the effect of PMF is not influenced by the time evolution of the cut 
off scale kc for this range of n, we approximately set E oc in the following analysis. 

III. CORRELATION IN POWER SPECTRA 

In this section, we define a power spectrum function of matter density with the PMF. In the linear approximation, 
solutions of eqs. (P5)l -(|33 p are divided into those with and without PMF by Green's function method as, 

(5(fc)- VL](fc)+ WlW- (34) 

Possible origins of PMF have been studied by many aurthors, however, we have not reached critical conclusions on 
the origin of PMF. Thus we cannot know how the PMF correlates with the primordial density fluctuations. However, 
almost all of previous works investigated the effects of PMF on density perturbations with the assumption that there 
is no correlation between them [23|. In order to study the PMF effect in a more general manner, we have to consider 
a correlation between the PMF and the primordial density fluctuations. Therefore we introduce "s" to parameterize 
the correlation between the PMF and the primordial density fluctuations. The power spectra of baryon {Ph{k)) and 
CDM (Pcdm(^)) density with the PMF in the linear approximation are. 



Pb(fc) = (^VFL](^)'^[b:FL](fc)/ + \VPMF](fc)'5*b:PMF](fc)/ 

+ 2^5[b:FL](fc)'^[b:PMF](fc)) ' (35) 

-PcDM(fc) = ^'^[CDM:FL](fc)'^['cDM:FL](^)) + ('^[CDMiPMF] (fc)<^*CDM:PMF] (^)) 

+ 2 {^(5[CDM:FL](^)'5[CDM:PMF](^)) ' (36) 

where, we define the cross correlations as, 

('5[b:FL](fc)'5*b:PMF](fc)) = ('^[b:FL](fc)'5*b^FL]('=)) ('^[biPMF] (fc)^*b:PMF] C^)) ' (3^) 

( '5[CDM:FL](fc)'5['cDM:PMF] 

(fc) ^ sM 

'^[CDMiFL] (^)'5[ci3M:FL] (^) ) ( '^[CDM:PMF] (^)'5[ci3M:PMF] 

(fc)), (38) 

where Jq, a G ([b : FL], [CDM : FL]) are baryon and CDM density fluctuations without the PMF respectively, and 
5(3, (i G ([b : PMF], [CDM : PMF]) are baryon and CDM density fluctuations with the PMF, respectively. When 
0<s<l,s = 0, and — 1 < s < on eqs. ([57)) and they stand for the positive, no, and negative correlations, 

respectively. 

The square root of power spectrum functions of Lorcnz force n[EM:S] C^: ''') in Eq. (|33[) does not have information about 
negative or positive, in other words, there is no information which of the magnetic pressure or tension is dominant, 
and whether the directions of forces from them are same or different. However, there must be such information, so we 
should take it into account. The Lorenz force term in Ea. p3p can be divided into two terms, the magnetic pressure 
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and the tension, of which ampUtudes are decided by eas. l|23|) and ([25]) . respectively. Comparing those equations, we 
can decide which of them is dominant in the Lorenz force term. To answer the first question we show in [1] which of 
the magnetic pressure or the tension dominates in the Lorentz force term. The former dominates when n < —1.5, the 
latter does when n > —1.5. 

As for the second question, we found that the roles of magnetic field pressure and tension in the Lorentz force 
term are different from each other for the random primordial magnetic field defined in Eq.(17). In other words, the 
scalar force from magnetic field tension acts on density field in the opposite direction in a statistical sense from what 
magnetic field pressure does. The reason is that, the variance of the Lorenz force |n(A:)| is always smaller than the 
simple sum of the variances of magnetic field pressure and tension, because the cross correlation between the two, 
^[EM:Si]2^[EM S2]' always gives the positive values for all k and Us considered in this paper. Since the force field from 
the magnetic field pressure is directly related to the magnetic field energy density distributions, that from the tension 
can also be related to the magnetic field energy distributions. Thus we can decompose the factor as 

S = S[LF] X Spp], (39) 

where 

"[LF] - I I, n> -1.5 (II), 

and 

< S[DF] < 1 (i), 
S[DF] = (ii), 
-1 < S[DF] < (iii). 

Here S[lf] represents (I) pressure dominant case, (II) tension dominant case. On the other hand, sppj represents 
(i) positive correlation between the matter and PMF distributions, (ii) no correlation, and (iii) negative correlation. 
Thus if s < 0, it means that the the matter and PMF distributions are positively correlated (spF] > 0) and the 
PMF pressure dominates in the Lorenz term (n < —1.5), or the matter and PMF distributions negatively correlate 
(■S[DF] < 0) and the PMF tension dominates in the Lorenz term (n > —1.5). In these cases the PMF effect acts against 
the gravitational collapse and makes the evolution of density purturbations slower like the gas pressure. On the other 
hand, if s > 0, the matter and PMF distributions positively correlate (spF] > 0) and the PMF tension dominates 
in the Lorenz term {n > —1.5), or the the matter and PMF distributions negatively correlate (s[df] < 0) and the 
PMF pressure dominates in the Lorenz term (n < —1.5). In these cases the Lorenz force of PMF accelerates the 
gravitational collapse. The above discussion is transparent especially in the perturbation evolution after decoupling, 
because 6 does not ocillate there. 



IV. NUMERICAL RESULTS 



In this section, we show our mumerical results. Since the recent upper limit on PMF amplitude by CMB is B\ ^ 7.7 
nG at ns = — 1.500l[l^. In order to be consistent with this result, we use the PMF parameter sets in the following 
calculations as, {Bx^ub) = (0.5 nG,-1.0), (0.5 nG, -2.001), (1 nG,-1.0), and (1 nG, -2.001) [4l|. 



A. Before decoupling 



The effect of PMF on CDM is much smaller than that on baryon before decoupling because the PMF cannot affect 
CDM directly. Moreover, the density of baryons oscillates with that of photons and their gravitational effects on 
CDM are very small. So we only consider the PMF effect on baryons here. 

Since the PMF can both effectively increase or decrease the gas pressure, the PMF affects the frequency of acoustic 
oscillations before the baryon fluid decouples from the photon fluid. This is important when considering the CMB 
anisotropics. Since the photon tightly couples with the baryon before the PLSS, it is natural that the frequency of 
oscillation of the photon density is indirectly increased or decreased by the presence of PMF. In this way the PMF 
affects the CMB photons which we can observe at present (see [TtI. ITsj). 
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B. After decoupling 

After baryons decoupling from photons, the baryons density evolution starts to affect CDM density evolution 
through gravitational interaction [3y|- So the PMF effect on the CDM increases with time through baryons (left 
panels in Figs. [Hand [3]). Once the CDM density fluctuations are generated by the PMF indirectly, they grow due to 
the gravitational instability so that the growthrate is the same as that of the primordial CDM density fluctuations. 
The baryons density fluctuations are generated by PMF directly, and therefore, the gravitational potential does not 

dominate the evolution of baryon density fluctuation until AirGp S ^ '^''^ "xpl^ '^"' • baryons density fluctuations 
with the PMF increase with the wavenumber k (right panels in Figs. [3] and [51). The reason is that the effect of the 
PMF is relatively larger at smaller scales (larger k values) in Eg. (p3|) . From Figs. |4]and[5]we found that the amplitude 
ratio between the density spectra with and without PMF {\P{k)/ Po{k)\ at k > 0.2 Mpc~^) lies between 75% and 
130% at present for the range of PMF parameters n = -2.001 and -1.0, 0.5 nG < i?;^ < 1.0 nG, and -1 < s < 1. 

V. DISCUSSIONS OF PMF MODEL 

Taking into account the stochastic PMF power spectrum sets (eqs.([20 |) -(f22 |) ). we can investigate the most accurate 
effect of PMF on the evolution of photons, baryons and CDM density fluctuation. From figsl2][5l it is turned out 
that there is a strong degeneracy of the correlation factor spp]; and the PMF amplitude B\. The correlation factor 
S[DF] depends on the PMF origin. So, in order to find a clue for solving such degeneracy problem, we shall discuss a 
relation of the correlation factor spp] and the PMF origin. 

A. Negative correlation of the density perturbations 

When spF] < 0, the tension of PMF delays the evolution of the matter density fluctuation, and the pressure of 
PMF accelerates it. In this case, therefore, when the PMF pressure for n < —1.5 dominates, the evolution of the 
matter density fluctuation is accelerated by the PMF. When n > —1.5, on the other hand, the evolution is delayed by 
the PMF. In order to substantiate such condition, the PMF must have been generated in the lower density regions. In 
order to substantiate such condition, the PMF must have been generated in the lower density regions. Some previous 
works mention that the PMF is directly proportional to the matter density on cosomological scale [s^ H^l- In these 
cases, since it is natural that there are more amplitude of PMF on the higher density regions in the frozen-in, this 
condision would be difficult to be realized with the causally generated PMF. 

B. No correlation 

In the case of the no correlation between the PMF and the matter density fluctuations, the power spectrum function 
of the matter denasity fluctuation with the PMF is increased by the PMF, independently dominances of PMF pressure 
or tension. Here we must notice the difference between the power spectrum function P{k) and the matter density 
fluctuation S. While the total density fluctuation S can be smaller or larger if the effect of PMF is dominated by its 
pressure or by its tension, respectively, the power spectrum function P{k) is always increased when PMF does not 
correlate with primordial density fluctuations. If the PMF is generated by density fluctuations proposed by [Sji peaks 
of PMF are at peaks of the pressure gradient of cosmic fluids. Consequently, the PMF is created along the border 
between high and low density regions, i.e., S ^ 0. In this case there would be no (or very weak) correlation between 
the PMF and density perturbation statistically. 

C. Positive correlation 

When S[DF] > 0, the pressure of PMF delays the evolution of the matter density fluctuation, and the tension of 
PMF accelerates it. So, in this case, when the PMF pressure dominates for n < —1.5 , the evolution of the matter 
density fluctuation is delayed by the PMF. In order to substantiate such condition, the PMF may be generated on 
the higher energy density regions. As an example, let us consider the PMF which was generated by a vector potential 
generated from the dilaton during inflationjj]. If the coupling between fields of dilaton and inflaton is negligibly 
small, there would be no correlation between the PMF and density perturbations. This is because the PMF was 
generated by the vector field coupled with dilaton, while the inflaton was responsible for density perturbations. 
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However, if we consider the case that a curvature coupling (Hke RF^^F^^^) generates the PMF in the same time, the 
positive correlation between PMF and density fields would be expected, since the electromagnetic fields are coupled 
with R or Hubble parameter, which is determined by the density fields. Because, as mentioned above Subsection lV A| 
the positive correlation between the magnetic field and the matter density at present is reasonable [s^ without 
any surprising PMF generations after the inflation, the positive correlation is the most natural consequence of such 
(inflationary) PMF generation models. 



VI. SUMMARY 

We numerically investigated the effect of the PMF on the energy density fields by considering the stochastic one that 
depends on scales, and we quantitatively discuss the effect of the PMF on the seeds of LSS in the early universe. We 
considered more general effects of the PMF than those considered in the previous works. We considered not only the 
magnetic field tension but also the increases of pressure and energy density perturbations from the field. Furthermore, 
by considering the correlation between the PMF and the matter density fiuctuation, and taking the mathematically 
exact stochastic PMF power spectrum sets, we obtained reasonable and accurate evolutions of baryon, CDM, photon, 
and therefore the large scale structure. We show that the PMF can play very different roles on the evolution of density 
perturbations accoding to the correlation. After decoupling, CDM is also influenced indirectly by the PMF through 
gravitational interaction. We have estimated the effects, and we found that the amplitude ratio between the density 
spectra with and without PMF (|P(A:)/Po(/c)| at /c > 0.2 Mpc~^) lies between 75% and 130% at present for the range 
of PMF parameters n = -2.001 and -1.0, 0.5 nG < Pa < 1.0 nG, and -1 < s < 1. 

Interestingly, it is reported that the magnetic field at large scales (A =lMpc) around B\ ^ nG [isl. [TtI. \\K Issj 
provides a new interpretation for the excess of CMB anisotropies at smaller angular scales. If the PMF with such 
strength was present, it is very likely that it has affected the formations at large scale structure as shown in the present 
paper. Yoshida, Sugiyama and Hcrnquist ^37| suggested that in order to avoid false coupling of the baryon and CDM 
for small scales, using independent transfer functions for the baryon and CDM is preferable. The PMF would be 
another source of this difference in the transfer function for baryon and CDM. Since the density perturbations in the 
early universe have evolved to the LSS at present age, the evolution of the LSS with the PMF becomes more different 
than that without the PMF. We have shown that the baryon and CDM energy density perturbations follow very 
different evolutions in the presence of the PMF; with the PMF taken into consideration, the evolution of large scale 
structure should become more complicated. 
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FIG. 1; Ratio of stocahstic PMF pressure and tension sources, y'(r[EM:Si] {t'^)T*^j^.g-i] (fe))/(r[EM:S2] (fc)T"[EM-s2] (^)) ' ^ function 

of wb- For illustration, a cut off scale kc is fixed to fcc = 10 Mpc~^. The pressure dominates for hb < —1.5, and the tension 
dominates for ub > —1.5. 



APPENDIX A: POWER SPECTRUM OF PMF 



In this section we derive power spectral of PMF, Z, 11 and E, which we used in sec|TIl The electro-magnetic 
energy momentum tensor can be decomposed into three parts, i.e., scalar, vector, and tensor parts. The scalar part 
of electro-magnetic energy momentum tensor Tij[EM:S] is definied as, 

Tij[EM:S] 'PijlmT''"\-EM]j (Al) 

where ^Vijim is a scalar project tensor: 
and 

py (fc') ^ (A3) 

Therefore it is easily shown that 

fc"' T'ijlm — kikfn- (-^4) 
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FJ =-2.tH]1 I 



Ji,-l.UaiG,ff= 1- 
J?i- I.OaCt, 



B^-(3.ijLQ.3- 1 ■ 
B,.-n.^iLfi, ■ 




Q.01 0.1 



0.01 0.1 



FIG. 2: Differences of transfer functions of CDM and baryon with (P(fc)) and witliout {Po{k)) the PMF, normalized by Po{k), 
with various strength of magnetic fields as functions of scale factor a at each wave number k as indicated. The curves in right 
and left panels correspond to the differences in CDM and baryon, respectively. The red thick, light-blue thin, orange normal, 
blue dotted bold, black dotted thin, and gree dotted middle curves in each both side panels show the differences for parameters 
(Ba,s) = (l.OnG,!), (1.0nG,0), (l.OnG,-!), (0.5nG,l), (0.5nG,0), and (0.5nG,-l), respectively. In aU figures, the power spectral 
index of the PMF is fixed to n = -2.001. The curves with s = in all panels are slightly above unity, although it is difficult to 
read off from the figure. 



1. Power spectrum of Lorenz Force: n(A;) 
Using Eqs. (|A2p - (|A4p . we obtain a two-point correlation function of scalar part 

(r[EM:S] ('s)^'[EM:S] (fc')) = k,k'kik^{Ti 

[EM]('s)T'm[EM](fc')) 

= (27r)3|n[EM:S](fc)l'5(fc-fc')- 



(A5) 



where we define a power spectrum of Lorenz force n[EM:S](^)- The electromagnetic stress-energy tonsor in k space is 
given by 



1 f (pk' n 



47ra4 J (2 



(A6) 
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FIG. 3: Same as Fig[2] but for a different power spectral index of the PMF n = -1.0. 
For convenience, we decompose 7j(fc)[EM] into two parts as follows, 

?7('^)[EM] = '7j(fe)[EM:l] - '7j(fe)[EM:2] 

1 r li^k' 



(A7) 
(A8) 

(A9) 



Correspondingly, we define TjemiSi] = Wk^ T- [em-.i] and T[em:S2] = fci^^^j [em:2]- Using Eq. (jA7l - lA9 p . we can rewrite 
the two point correlation function of Lorenz force of scalar part as, 



(7'(fc)[EM:S]T*(fc')[EM:S]) ((r(A;)[EM:Sl] " T^l^) [EM:S2] ) (7"* (fc') [eM:S1] " r*(fc')[EM:S2])) 

= (r(A;)[EM:Sl]7'*(fc')[EM:Sl]) " (r(fc)[EM:Sl]r*(fc')[EM:S2]) 
- (r(fc)[EM:S2]r*(fc')[EM:Sl]) + (r(fc)[EM:S2]r*(fc')[EM:S2])- 



(AlO) 
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n = -2.m] I 





i?i-4J.5aG, J- 


Bi.= LtiiiG,i= 




Bv - \.(SnQ,fi-~l 


ft: - O.^fiG, A - - 1 




jtij' iiiri fur H. ■: f|.*pO 




.0.05 O.J St. 1 5 0.25 0,3 



0.1:5 



FIG. 4: Differences of transfer functions of CDM and baryon with (P(fc)) and without (Po(fe)) the PMF, normalized by Po{k) 
as a function of wave number k at each redshift as indicated. The curves in right and left panels correspond to the differences 
in CDM and baryon, respectively. The red thick, light-blue thin, orange normal, blue dotted bold, black dotted thin, and 
gree dotted middle curves in each both side panels show {B\,s) — (l.OnG,!), (1.0nG,0), (l.OnG,-!), (0.5nG,l), (0.5nG,0), and 
(0.5nG,-l), respectively. In all figures, the power spectral index of the PMF is fixed to n = -2.001. The curves s = in all 
panels look like unit, acutually, those is more than unit. 



Also, from (IA5[) . n[EM:S]('^) Eg. ([^5)) and ([55]) can be written by, 

n[EM:S](fc)(5(fc- fe') = W^^(T(fc)[EM:S]r*(fc')[EM:S])- 



(All) 
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FIG. 5; Same as FiglH but for a different power spectral index of the PMF ns ~ -1.0. 



2. Power spectrum of Shear Stress: Z{k) 

A power spectrum of Loreze force for the scalar part is decomposed into trace-trace and traceless (a shear pressure 
of PMF) parts. A traceless part of electromagnetic stress-energy tensor can be written by, 



-J [EM] 
1 



t: 



J [EM] 



Ana'^ J (27r) 



g'^J^fc [EM] 



-5]B\k')Bi{k - k') - B\k')Bj{k - k') 



(A12) 
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Thus, a scalar part of traceless component is given by, 

1 f cPk' f - fl 



47ra4 7 (27r) 

- |^r(fc)[EM:Sl] - T(fc)[EM:S2]| (A13) 

Similarly as the above subsection, we obtain the two point correlation function of PMF shear stress for scalar part as 
the following, 

(S(fc)[EM:S]S(fc')[EM:S]) 

4 2 

= g('7"(^)[EM:Sl]^*('^ )[EM:S1]) ~ g (^(^) [EM:S1]^"* ) [EM:S2] ) 

2 

-2(^'(fc)[EM:S2]T'*(fc)[EM:Sl]) + (r(A;)[EM:S2]7'*(^)[EM:S2]) 

= {2nf\Z{k\^M:S]?5{k-U), (A14) 
where 2'(fc)[EM:S] is a power spectrum of PMF shear stress. We rewrite Ea. (|A14p as following, 

Z2(fc)[EM:S]'5(fc- fc') = ^(S(fc)[EM:S]S(fc')[EM:S])- CA15) 

3. Power spectrum of PMF energy:i5(fc) 
The PMF energy is defined by Eq.(|3]). The PMF energy in k space is given by 

ro°(fc)[EM] = j2^B\U)Bi{k-U) = r(fc)[EM:Sl]. (A16) 

Similarly as the above subsection, we obtain the two point correlation function of PMF energy for scalar part as the 
following, 

(r(fc)[EM:Sl]r(fc')[EM:Sl]) = {'2^f\E{k\^M:S]?5{k-U), (A17) 

where i?(fc)[EM:S] is a power spectrum of PMF energy. We rewrite Eq. (jA17p as following, 

E\k\^M:S]5{k-U) = ^(r(fc)[EM:Si]r(fc')[EM:Si])- (A18) 

In order to obtain n(fc)[EM:S], 2'(fc)[EM:S], and i^(fc)[EM:S], we calculate (r[EM:Si] (fc)T[EM-si] (^)) > 

(r[EM:Si](A;)rfEM:S2](^)) + (T[EM:S2](fc)TfEM:Si](^)) ' and (T[EM:S2](fc)T*j,M:S2](fe)) i^ Eqs-M), jMl and dSni the 
following subsections. 

4. {7'[EM:Sl](fe)T[EM:Sl](fe)) 

Using Eg. ljAip . (jA2[) and (|A8p . the scalar part of T(fc)[EM:i] becomes 

T'(fc)[EM:Sl] = kk^T^{k)[EM:l] 



(A19) 
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Thus the two point correlation function of Eq. (|A19[) given by 

(T(fc)[EM:Sl]7'(p)[EM:Sl]) = 

^4(2^)8^8 / d'k'd'p'{B\k')B,{k~ k')B'*{p')B:{p- p')) 

We assume the random magnetic field is Gaussian and apply the Wick's theorem 

(^B\k')B,{k-k')B'*{p')Bi*{p-p')) = {B\U)B,{k-U))l^B'*{p')Bi*{p-p') 

+ (^B%k')B'*ip')) {B,ik - k')Bi*ip - p')) 
+ {B\U)Bi*{p - p')) (^B,{k-U)B'*{p') 

A two-point correlation funciton for PMF is defined by Ref.(l3| 



{B\k)B,{k-k!))(^B' {p')Bap-p')) = y-:^Y{^\ k"\~prm)5{-p), 



and 



B\U)B'\p')){B,{k-U)Bi*{p~p')) = k'-\k-U\^^P^\k')Pa{k-k') 

xS{k' -p')d{{k-k')-ip-p')), 



{B'\U)Bnp-p'))(B,{k^U)B'*(p')) = k-^\k-U\-Pl(k')Pl(k-k') 



x6{{k-k')-p')Sik' -ip-p')), 



(A20) 



(A21) 



B\k)B=\U)) = ^-±__^k"P^^ik)Sik-k'). (A22) 
Using the reality condition 

B\k)^ B'*{-k), (A23) 
and substituting (jA22| into (|A2ip . first, second, and thrid terms become 

/■ \ 2 

^r,\_)y^T^) I -'^^\5{k)5{-p), 

(A24) 



(A25) 



(A26) 



respectively. Because k ^ Q{x = 27r/fc ^ oo), 5{k) = S{—p) = 0, the first term on the r.h.s.in Eq. (jA21[ ) is 0. From 
Eq.(jA3| 

P''(mKfc-fcO = l+ ^'^,,(^':'y . (A27) 
Therefore second and third terms in Eq. (jA21|l become 
'B\U)B'*{p')) {B,{k-U)Bi*{p-p')) 

(A28) 
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and 



1 J 


\ (27r)"+s 


Bl 


23(27r)7a8 ] 


[ 2fc^'+^ r 


(^) 




^^y' dc\k 











(A29) 



(A30) 



{B\h!)Bi*{p-p'))I^B,{k-U)B'\p')) 
respectively. Using Eq. (|A28|) and (|A29p . Eq. (|A30p becomes 

(T(fc,T)[EM:Sl]7'*(P, 'r)[EM:Sl]) 

X {(5(fc' - p')5{{k - k') -{p- p')) + 5{k! -{p- p'))5{{k - k') - p')} 
Integrating Eq. (jA30p by p' , we obtain following equation 

{T{k,T)[EM:Sl]T{p, r)*EM:Sl]) 

We define that 

C = cose = fe- = -^-7^. (A32) 

fc'K 

Substituting Eq. (jA32| into Eq. (|A3ip , The two point function of Lorenz Force becomes 

1 f (27r)"+* 1 ^ 

(T(fc,T)[EM:Si]r(p,r)['EM:Si]) = 23(2T)8n8 I 2/;;+3 Fp^J 

X y" (fk'k"'\k- Uy^-"^ {(1 + C2)/c2 _ Akk'C + 2k'^} S{k - p). (A33) 

Choosing k to the polar axis as 

Sk' = k''^dk' sin c dc d(t>, (A34) 
and integrating Eq. (|A33p by (p, the two point function of Lorenz Force is given by 

{T{k, t)iem:Si]T{p, t)[em:Si]) 



(A31) 



(A35) 

Almost all of the previous works have treated terms which include C in the middle parenthesis as unity. In this paper, 
however, we calculate Ea. (|A35p furthur by integrating by parts. In addition they have calculated integration of k' 
using the Taylor expansion by k' /k(k' <C k) or k/k'{k ^ k'). If we would want to estimate Eq. (jA35[) for only k' <^k 
or k <^ k' , such approximation would be useful. However there is the value A:' ^ fc in the integration range, so we 
must caluculate Eq. (jA35|) without such Taylor expansion. Integrating Eq. (jA33[) by C, after long but straightforward 
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calculation, we can obtain following expression, 

1 f (27r)"+^ B'^ Y 
(r(fe,r)[EM:Si]r(p,T)[EM:Si]) = 4(2T)7a8 I 21:'^+^ T ("+^)" J 

+ 4n + 1 



1 



fcfc'n(n + 2)(n + 4) 
{|fc-A:'|"+2 + |A: + fc'|"+'} 



{(fc + fc')"+'-|fc-fc'|"+'} 



k'^n{n + 4) 



+ TTw ^^T7 JT {(k + k'r^' - \k - 



(A36) 



5. (7'[EM:Sl]r[EM:S2]) + (2"[EM:S2] TfEMiSl] ) 

Using Eas. (|Al[) . \k2\ and (|A9|1 . the scalar part of T(fc)[EM:2] becomes 

7'(fc,T)[EM:S2] = r)[EM:2] 

Thus cross correlations between Eqs. (|A19p and (|A37p are given by 

(T'[EM:Sl]('s)'7'[EM:S2](P)) = 

j d'k' J d'p'kik"\B^k')B,ik~k')B'*ip')B*^ip-p')), 



23(27r)8a8 
and 



(ASS) 



(%M:S2](fc)7'[*EM:Sl](P)) 
1 



23(27r)8a8 

respectively. Using Eq. (IA2ip . Eq. (|A38p and (|A39|) become 

1 f(27r)"+8 B2 ^2 



d^k' / d^p'hkP{B'{k')Bj{k-k')B'i*{p')B*{p~p')), 



(A39) 



(T[EM:Sl](fc)TfEM:S2](P)) - ^3(^ j / / '^'^'j ^'p' k'-\k - U^p^p"^ 

X {P«'(fc')P,™(fc - - p')5{{k - k') -{p- p')) + P\{k- k')PUk')S{k' -{p- p'))S{{k - k') - p')} 

(A40) 



and 

\n+S r2 ^ 



(r[EM:S2](fc,r)r[*EM:si](P,^)) = 23(2^)8a« { ^2fcf+^ r(^) | / ^'^'/ rfVfc'"|fc- fcrfc.fc^P™(fc')P.,(fc- fc') X 
{5(fc' - p')'5((fc - k') -{p- p')) + S{k' -{p- p'))d{{k - k') -p')}, (A41) 
respectively. Integrating Eq. (jA40p and Ea. (|A41[ ) by p', we obtain following equations, 

{T[EM:Sl]{k)T^EM:S2]iP)) = 

1 f ('27r'l"+^ 1 ^ /■ 

(A42) 
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and 

('7'[EM:Sl]('i:)?'[EM:S2]('^)) = 

1 f(27r)"+8 B 



, 2 

2 



(A43) 



respectively. Similary as above the subsection, term containing the unit vector k in Eas. (jA42p and (jA43p becomes 
Substituting Ea. (|A44p into Eqs. (|A42p and (|A43|) .and combining Eqs. (|A42|) and (|A43|) . we obtain following equation 

y" d=^fc'fc'"ifc- fc'i"-2^-'(i - c^){k' - kc)s{k-p). 

(A45) 

Similary, integrating Eq. (jA33P by C, after long but straightforward calculation, we can obtain following equation, 

' (27r)"+8 Bj 



1 j (27r)"+8 Bl 



2 -2 



1 (27r)"+'* i?^ /■ , 

(r[EM:Sl](fc,T)r[EM:S2](fc.^)) + (%M:S2](fc,T)TfEM:Sl](fe^T)) ^ (27r)7a8 | 2/;"+^ r("+^) J '^'^^ ~ J 

^ {(fc + - (fc - - ,,.,3 , — -- {(fc - fc')"+* + {k + kT+^} 



I 1/71 + 3 



{kk'fn{n + 2) ' V ; J pfc'3^(„ + 2)(n + 4) 

^ {(fc + fc')"+' - (fc - A;r+'} +TT77r7 ^1 Wi ^ {(^ + ^')"+' - (fc - 



A;3A:'2n(n + 2)(n + 4) ' ' 'J A:3fc'4n(n + 2)(n + 4)(n + 6) 

(A46) 

6. (r[EM:S2]r[EM:S2]) 

Using Eq. (jA37l )■ a two point correlation function of T\^EM:S2]{k) becomes 

(r[EM:S2](fc)rfEM:S2](p)) = ^2(^^ / k' p' W Pi^ {B'{k')B,{k - k')B'*{p')B„:{p-p')) (A47) 
Useing Eas. (jA2ip and (jA22p .forthremorc integling by p\ we get 

(r[EM:S2] (fc, t)7'*eM:S2] (P. "t)) 

= ^(^'P) 22(27r)«a« | 2^+^ T {^) \ J d'k'k'"\k~ k'\"hypir{P^\k')P^"^ik - k') + P.Uk')P',{k - ^0} 

(A48) 

Similary as the above subsection, the term containing the unit vector k on Eq. (jA48|) is expressed as 

kk'pif''P'\k')P,^(k - k') = (1 - C')' . (A49) 

Similary, substituting Ea. (|A49[) for Ea. (|A48p , we obtain following equation 

(%M:S2] {k-. t)7'[*EM:S2] C^' ^)) 



2(27rra8| 2fc^3 T (4^) } J ^k' J_^dCk'-+^\k k'\ 



(A50) 
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Similary, integrating Eq. (jA50p by C, after long calculating, we can obtain the following expression, 

1 f (27r)"+^ \^ f I I 4 r / / 

"(2^0^ I 2fc;:+3 r(4^)j y '^^'^'"^'(fcfc')'«(n + 2)(n + 4) i {(^ + fc')"+'-|^-^r+'} 

(A51) 



APPENDIX B: CUT OFF SCALE OF PMF 



The photon mean free path is 

aTne[T) 

where rie is the electron number density. From Refs. [H, [2^ , the cut off scale of PMF is defined by 



k^^(T) = <j (o^ ' ^ ^ (B2) 

[ k(j (Tdcc), T > Tdoc 

where Tdcc is the time of decoupling of photons from baryons. Following Ref . [l3| . for a stchastic magnetic field with 
a power-law power spectrum, the relation between the effective homogeneous field B and B\ is written by 

rt + 3 

kc ^ 



B = B^i^^] . (B3) 

Subtstituting Eq.(|B2| into Eq.(|B3|, we obtain 



fcc W = i 47r(p+p| Jo T' ' ^ 'doc ^g4) 
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